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ABSTRACT 

Context. Hot subdwarf-B stars in long-period binaries are found to be on eccentric orbits, even though current binary-evolution theory 
predicts those objects to be circularised before the onset of Roche-lobe overflow (RLOF). 

Aims. We aim to find binary-evolution mechanisms that can explain these eccentric long-period orbits, and reproduce the currently 
observed period-eccentricity diagram. 

Methods. Three different processes are considered; tidally-enhanced wind mass-loss, phase-dependent RLOF on eccentric orbits 
and the interaction between a circumbinary disk and the binary. The binary module of the stellar-evolution code MESA (Modules 
for Experiments in Stellar Astrophysics) is extended to include the eccentricity-pumping processes. The effects of different input 
parameters on the final period and eccentricity of a binary-evolution model are tested with MESA. 

Results. The end products of models with only tidally-enhanced wind mass-loss can indeed be eccentric, but these models need to 
lose too much mass, and invariably end up with a helium white dwarf that is too light to ignite helium. Within the tested parameter 
space, no sdBs in eccentric systems are formed. Phase-dependent RLOE can reintroduce eccentricity during RLOF, and could help 
to populate the short-period part of the period-eccentricity diagram. When phase-dependent RLOF is combined with eccentricity 
pumping via a circumbinary disk, the higher eccentricities can be reached as well. A remaining problem is that these models favour a 
distribution of higher eccentricities at lower periods, while the observed systems show the opposite. 

Conclusions. The models presented here are potentially capable of explaining the period-eccentricity distribution of long-period sdB 
binaries, but further theoretical work on the physical mechanisms is necessary. 
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1. Introduction 


Hot subdwarf-B (sdB) stars are core-helium-burning stars with 
a very thin hydrogen envelope (Mh < 0.02 Mq), and a mass 
close to the core-helium-flash mass ~ 0.47 Mq ( |SafFer et al.| 
1994[ [Brassard et al.|[200T| ). These hot subdwarfs are found in 
all galactic populations, and they are the main source for the UV- 
upturn in early-type galaxies (Green et al.|1986||Greggio & Ren-j 


zini||1990] [Brown et al.||1997 1 . Furthermore, their photospheric 


(|Webbink|1984|l.The possibility of forming a single sdB after a 


chemical composition is governed by diffusion processes caus¬ 
ing strong He-depletion and other chemical peculiarities ( [Heber[ 
1998| l. The formation of these extreme-horizontal-branch objects 
is still puzzling. To form an sdB star, its progenitor needs to lose 
its hydrogen envelope almost completely before reaching the tip 
of the red-giant branch (RGB), so that the core ignites while the 
remaining hydrogen envelope is not massive enough to sustain 
hydrogen-shell burning. A variety of possible formation chan¬ 
nels have been proposed. Currently, there is a consensus that sdB 
stars are only formed in binaries. Several evolutionary channels 
have been proposed, where binary-interaction physics plays a 
major role. Close binary systems can be formed in a common- 
envelope (CE) ejection channel ( Paczynski||1976 1 , while stable 
Roche-lobe overflow (RLOF) can produce wide sdB binaries 
( [Han et al.|2000||2002[ l. An alternative formation channel form¬ 
ing a single sdB star is the double white-dwarf (WD) merger, 
where a pair of white dwarfs spiral in to form a single sdB star 


merger with a sub-stellar object has also been proposed (Soker 
2014| l, but not modelled. This channel could produce sdBs with 
a narrow mass distribution unlike the WD mergers. 

Han et al. (2002 2003| l addressed these three binary- 
formation mechanisms, and performed binary-population- 
synthesis (BPS) studies for two kinds of CE-ejection channels, 
two possible stable-RLOF channels and the WD-merger chan¬ 
nel. The C- ejection channels produce close binaries with pe¬ 
riods of Porb = 0.1 - 10 d, and main-sequence (MS) or white- 
dwarf (WD) companions. The sdB binaries formed through sta¬ 
ble RLOF have orbital periods ranging from 10 to 500 days, and 
MS companions. Chen et al. (|2013 i revisited the RLOF mod¬ 


els of Han et al. (2003[) wTS more sophisticated treatment of 


angular-momentum loss. When including atmospheric RLOF, 
these revised models can reach orbital periods as long as ~1600 
d. Finally, The WD-merger channel can lead to single sdB stars 
with a higher mass, up to 0.65 Mq. A detailed review of hot sub¬ 


dwarf stars is given by Heber (2009 1 . 


Many observational studies have focused on short-period 


sdB binaries (Koen et al.|1998[|Maxted et al.|20001|2001[|Heber| 

et al.|2002| Morales-Rueda et al.|20031 

Napiwotzki et al.|2004 

Copperwheat et al. 20111, and over 100 

3f these systems are cur- 
dix A). These observed 

rently known (jGeier et al.|2011 Appen 


short-period sdB binaries correspond very well with the results 
of BPS studies. Currently, only nine long-period sdB binaries 
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are known ( 

Green et al. 

2001 

0stensen & Van Winckelj 

2011 

2012 

|Deca et al.||2012 

IBarlow et al.||20121 |Vos et al.| 

2012 

2013 

[20141. Even though this is a small sample, their period- 


eccentricity distribution is not compliant to the prediction of evo¬ 
lution models. All current models predict circular orbits, while 
seven out of nine observed systems have a significantly eccentric 
orbit. 

This eccentricity problem is not entirely new, and a few pos¬ 
sible solutions have been proposed, although none have been ap¬ 
plied to the case of sdB binaries. There are three potential mech¬ 
anisms that can create eccentric orbits described in the literature. 
|Soker| ( [2000| l based on theoretical work of [Eggleton] ( |2006| pro¬ 
posed the mechanism of tidally-enhanced wind mass-loss to al¬ 
low the eccentricity of the orbit to increase. In this framework, 
the wind mass-loss is increased due to the tidal influence of a 
companion in an eccentric orbit, and the effect is present be¬ 
fore the system comes into contact. It reduces the tidal forces 
by keeping the sdB progenitor within its Roche lobe, while the 
phase-dependent wind mass-loss can increase the eccentricity. 
This mechanism is also known as Companion Reinforced Attri¬ 
tion Process (CRAP). |Siess et ak] P014| l used this mechanism 
to successfully explain the orbit of the highly eccentric He-WD 
binary IP Eri. 

A second possible mechanism is that of phase-dependent 
mass loss during RLOE Similar to the phase-dependent mass 
loss in the CRAP mechanism, a varying mass-loss rate in a bi¬ 
nary that is larger during periastron than apastron, can increase 
the eccentricity of the orbit. The difference with the previous 
mechanism is that it is active during RLOE. This mechanism 
was used by |Bonacic Marinovic et al.] ( |2008[ l, who used a model 
with enhanced mass loss from the AGB star due to tidal inter¬ 
actions and a smooth transition between the wind mass loss and 
the RLOE mass loss regimes to explain the eccentric population 
of post-AGB binaries. 

The third method of increasing the eccentricity of binaries is 
by interaction of the binary with a circumbinary (CB) disk. The 
motivation to include this process in the context of sdB stars, 
is the observational finding that stable circumbinary disks are 
commonly observed around evolved post-AGB binaries (e.g. de 
|Ruyter et al.||2006l |Hillen et al.||2014| and references therein). 
The longevity of these disks are corroborated by the strong pro¬ 
cessing of the dust grains as attested by the infrared spectral dust- 
emission features (e.g. |Gielen et al.|2008||201 1) and the millime¬ 
tre continuum fluxes that indicate the presence of large grains ( [de] 
Ruyter et al.|2005|l, while interferometric techniques are needed 


to resolve them (see e.g. [Hillen et al.||2014| and references 
therein). The Keplerian rotation is, up to now, only spatially re¬ 
solved in two objects ( jBujarrabal et al.|2005 2015|l. Recent sur- 
veys of the Large and Small Magellanic Clouds ([van Aarle et al.| 
[2011 1 [Kamath et al.||2014| ) show that a significant fraction of 
post-AGB stars show the distinctive near-IR excess indicative of 
a stable disk. One of the results of [Kamath et al.| ( |2014||20T5] l 
is that a significant population of post-RGB stars were identified 
with circumstellar dust likely in a disk as well. The disk evo¬ 
lution will determine the infrared life time of the systems and 
hence the detectability. While there is ample observational ev¬ 
idence that disks are common in evolved binaries, we assume 
here that disks were also present during the RGB evolution of the 
sdB progenitor. [Dermine et al.[ ( [2013[ l explored the eccentricity¬ 
pumping mechanisms of CB disks in post-AGB binaries, based 
on theoretical work and smooth-particle-hydrodynamic (SPH) 
simulatio ns of|Artymowicz & Lubow ( 1994| l and Lubow & Arty- 
[mowicz[ ( [1996[ ). An issue with the results of Dermine et al. 
(2013|l is that there was a mismatch between the disk mass dis- 
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Table 1. The observed periods and eccentricities of all known long- 
period sdB binaries. See also Fig.[T] 


Object name Period (days) Eccentricity 


PG 1701+359'’ 

734 

± 

5 

0.00 

± 

0.04 

PG 1104+243'’ 

753 

± 

3 

0.00 

± 

0.01 

PG 1018+243" 

759 

± 

3 

0.05 

± 

0.01 

PG 1449+653" 

909 

± 

5 

0.11 

± 

0.04 

Eeige 87'' 

936 

± 

2 

0.11 

± 

0.01 

BD+34°1543'' 

972 

± 

2 

0.16 

± 

0.01 

BD+29°3070'' 

1283 

± 

63 

0.15 

± 

0.01 

BD-7°5977" 

1314 

± 

8 

0.15 

± 

0.01 

Bal 8280003" 

1363 

± 

25 

0.25 

± 

0.02 


Notes. <° MBarlow et al.U2012b ^^MVos et al.U2012t , [Peca & Vos[ 


i 2015[ l , <‘') |Vos etal.| ( |26i3t ,*"> 


Vos et al. 


(20141 


tributions used to derive different parts of the interaction mecha¬ 
nisms. 

In this article we will explore these three methods to find for¬ 
mation channels that can explain the eccentricity of sdB binaries, 
using the stellar/binary evolution code MESA. In the case of CB 
disks that are formed during RLOE, the eccentricity pumping ef¬ 
fect of the disk is combined with the effect of phase-dependent 
RLOE. We do not aim to perform a binary-population-synthesis 
study, but to explore in a limited number of initial conditions, 
the efficiency of the different processes. Our aim is to describe 
the effect of the model parameters on the final period and eccen¬ 
tricity of the binary, and to discover which areas in the period- 
eccentricity diagram can be covered by the three methods. 

An overview of the currently known systems is given in Sect. 
[2 The MESA evolution code is explained in Sect. [^ The mod¬ 
elling methodology is explained in Sect.ffl While the different 
models are presented in Sects. [^ [^ and [T] for respectively the 
CRAP models, phase-dependent RLOE models and models con¬ 
taining a CB disk. The obtained period-eccentricity distribution 
is discussed in Sect.[^ and a summary and conclusion is given in 
Sect.[^ A detailed overview of the binary physics used in MESA 
is given in appendix [A[ 


2. Observed period-eccentricity diagram 


There are currently two wide sdB binaries with circular or¬ 
bits and seven significantly eccentric systems covered in the 
literature. The first observed long-period sdB binary was 
PG 1018-047 with a period of 760 + 6 days and an unclear ec¬ 
centricity ( Deca et aL[[2012| l. Based on new observations Deca 
[& Vos[ ( [2015[ l revised this to a period of 759 + 3 days and an 
eccentricity of 0.05 + 0.01. [Barlow et alT[ ( [2012[ ) published the 
orbits of two sdB systems, PG 1701-1-359 and PG 1449-1-653, the 
first circular and the second with an eccentricity of 0.11+0.04. 
There is one extra circular system known, PG 1104+243 ([Vos| 
et al.|2012[ ). Eurthermore, five eccentric systems were observed 
with the Mercator telescope (Vos et al. 2013| [2014[ ) using the 
HERMES spectrograph ( [Rasldn et al.[[20lT l. All known orbits 
are given in Table [T] while the resulting period-eccentricity di¬ 
agram is shown in Eig. [^ Even though it is based on published 
data, this distribution has not been shown this clearly before. 

Of the systems with the shortest periods, two have a circu¬ 
lar orbit (PG 1104+243, PG 1701+359) and one has a slightly 
eccentric orbit (PG 1018+243). Models should thus account for 
both circular and eccentric systems. Eurthermore, there is a clear 
trend detected in Eig.[2in which the eccentricity increases with 
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Period (days) 

Fig. 1. The observed period-eccentricity diagram of all known long- 
period sdB binaries. There is a clear trend visible of higher eccentricities 
at higher orbital periods. See also Table[T] 


increasing orbital period, ranging from e - 0.05 at P = ~750 d 
to e = 0.25 at P = ~ 1350 d. 


3. MESA 

Modules for Experiments in Stellar Astrophysics (MESA^is 
an open-source state-of-the-art ID stellar-evolution code, which 
amongst others, includes a binary module to compute evolu¬ 
tionary tracks of binary stars. The stellar-evolution modules of 
MESA are extensively described in the two instrument papers; 
[Paxton et al.| ( |20lT||2013| l. In this article we have used version 
7211 of MESA. A reason to work with this evolution code is the 
availability as open source. Moreover, the code calculates stellar 
models of low-mass stars through the helium flash. 

The binary module of MESA is under continuous develop¬ 
ment. We started from version 7211 and extended it with several 
physical processes necessary for this research. The most impor¬ 
tant of which are: 


Accretion of the mass lost in stellar winds. To calculate the 
accretion rates, the Bondy-Hoyle formalism as described in 
Hurley et al. ( 2002| l is used. 

Tidally-enhanced mass loss. Wind mass loss can be en¬ 
hanced by the tidal influence of the companion star. The 
Companion Reinforced Attrition Processes (CRAP) mech¬ 
anism of jTout & Eggleton] ( |1988| l is used to calculate this. 
Phase-dependent mass loss during RLOE or phase- 
dependent mass loss through stellar winds can possibly in¬ 
crease the eccentricity of the orbit. The formalism of jSokerj 
( |2000| l and [Eggletonj ( |200^ is used to determine the change 
in eccentricity due to mass lost to infinity and mass accreted 
by the companion. 

Circumbinary (CB) disks. Due to Lindblad resonances, CB 
disk-binary interactions can change the orbital period and ec- 

( fT994l l 
ate the 


centricity. The formalism of Artymowicz & Lubow 
and Lubow & Artymowicz ( 1996| l is used to calcu. 
CB disk-binary interaction. 


* MESA is available on: http: //mesa. sourceforge. net/ 


A list of the processes implemented in the binary module is 
given in appendix We limited ourselves to listing the phys¬ 
ical processes. Eor the exact software implementation, we refer 
the reader to the open-source code, and future instrument papers. 
Most of the binary additions summarised above, will be added to 
MESA in the future. Until then, a copy of the binary module used 
here will be made availabl^B 

As an example of binary tracks leading to post-RGB evo¬ 
lution, the evolutionary tracks in the HR diagram for a late hot 
flasher and an early hot flasher are displayed in panels B and C 
ofEig.g 

In this article we will refer to the sdB progenitor as the mass¬ 
losing star or donor star, using the subscript d. The companion is 
also called the accretor, even when it won’t accrete mass, using 
subscript a. 


3.1. Stellar input parameters 


As the main focus of this contribution is the study of binary- 
evolution processes, we will use standard parameters for the 
evolution of the individual stellar components. These include 
a standard atmospheric boundary condition at an optical depth 
of T = 2/3, a mixing-length parameter of auLT = 2 and 
default opacity tables (OPAL type i). For the sdB progeni¬ 
tor we used an extended version of the standard nuclear net¬ 
works to include all reactions for hydrogen and helium burn¬ 
ing (pp_cno_extras_ol8_ne22 .net). The initial composition 
is X = 0.68, Y = 0.30 and Z = 0.02. Furthermore we used a 
Reimers wind on the RGB with /yreimeis = 0-7 ( |Reimers||1975| l, 
and a Blocker wind scheme on the post-EHB with r/biocker = 0-5 
( |Blocker||l995| l. All models are calculated from the pre-main 
sequence till the WD-cooling curve, with the RGB taking the 
longest computation time. The sdBs resulting from these stan¬ 
dard settings do not completely reproduce the observed effective 
temperature or surface gravity . This has been reported before, 
for example by 0stensen et al. ( |2012 l; 0stensen ( [20141 ). We are 
aware of the discrepancy, but solving this is beyond the scope of 
this article. 


4. Modelling methodology 

The focus of this article is the effect of the different eccentricity¬ 
pumping mechanisms on the evolution. Three methods are inves¬ 
tigated in the following sections. Tidally-enhanced wind mass- 
loss is considered separately from the other two mechanisms. 
The main argument is that tidally-enhance wind mass-loss in¬ 
creases the wind mass loss so that the donor star never Alls 
its Roche lobe. In Sect.|^ the effects of firstly phase-dependent 
RLOE alone, and secondly phase-dependent RLOE in combina¬ 
tion with a CB disk are described (Sect.|7]i. As the CB disk is 
created from mass lost during RLOE, the eccentricity-pumping 
effects of phase-dependent RLOE will also be present when 
the CB disk-binary interaction starts. We do not consider the 
CB disk-binary interactions separately, but assume the disk is 
formed by mass lost in the previous phase. For the phase- 
dependent RLOE models and the CB-disk models, a default 
model is used to show the evolution of several binary param¬ 
eters. After which the effect of both initial binary parameters 
and method dependent parameters are discussed. The complete 
range in period an eccentricity that can be covered using any of 
the eccentricity-pumping mechanisms is discussed separately. 

^ The binary module is available on: http://ster.kuleuven.be/ 
~j orisv/MESA_binary 
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Fig. 2. HR diagram of the evolution of the sdB progenitor for different values of the initial period. The location of the He ignition is indicated with 
a red square. Panel A: P = 600 d, the donor star loses too much mass to ignite He and ends up as a cooling white dwarf (M^ = 0.448 Mq). Panel 
B: P = 650 d, a late He flasher (M^ = 0.456 Mg). Panel C: P = 750 d, an early He flasher (M^ = 0.466 Mq). Panel D: P = 800 d. The core is too 
massive, and the donor ignites He on the tip of the RGB (M^ = 0.929 Mq). See section|^for discussion. 


In the case of phase-dependent RLOF, and the CB disk- 
binary interaction models, tidal forces will have circularised the 
orbit before the onset of RLOF. However, if the eccentricity is 
zero, the eccentricity-pumping effect of phase-dependent mass 
loss and CB disk-binary interactions will be zero as well. There¬ 
fore, we impose a lower limit on the eccentricity of 0.001. This 
lower limit is chosen small enough to be reasonable. This min¬ 
imum eccentricity is applied in all models, also those with only 
tidally-enhanced wind mass-loss. Thus if we refer to circularised 
models, this effectively means models that reached the minimum 
eccentricity of e = 0.001. 


For each eccentricity-pumping method a substantial number 
of models with varying initial parameters and process-depending 
parameters are calculated. The process-dependent parameters 
and their ranges are discussed in the separate sections. The initial 
binary parameters (period, eccentricity and mass of both compo¬ 
nents) are determined based on the type of systems that we want 
to produce. We are focusing on the long period sdB binaries, 
thus initial orbital periods vary between 500 and 900 days. The 
sdB progenitor mass varies between 1.0 and 1.5 Mq. While the 
initial companion mass varies between 0.8 and 1.45 Mq, and is 
always lower than the sdB-progenitor mass. Apart for the tidally- 
enhanced wind mass-loss models, all models circularise com¬ 
pletely before the eccentricity-pumping mechanisms become ef¬ 
fective. Changing the initial eccentricity for these models then 
only affects the initial orbital momentum of the binary, thus 
changing e, has a similar effect as changing the initial period. 


Only a small subsection of this parameter space will result 
in binaries containing an sdB component. In the discussion of 
the effect of the process-depending parameters, only the models 
containing an sdB component are used. The selection criteria for 
sdBs are based on the stellar mass at He ignition, which has to 
be below 0.55 Mq. And the absence of a hydrogen burning shell. 
In our models, there is a clear mass gap between the sdB models, 
and the models that ignite He on the RGB. The latter having final 
masses of Md > 0.7 Mq, while the sdB models have final masses 
of Md < 0.49 Mq. Models that don’t ignite He are obviously not 
sdBs either. 


Article number, page 4 of 


18 


5. Tidally-enhanced wind mass-loss 


|Siess et aL] ( |2014[ ) has modelled the long-period eccentric system 
IP Eri by using tidally-enhanced wind mass-loss in combination 
with the eccentricity-pumping effect of phase-dependent mass 
loss. By increasing the wind mass loss, the radius of the donor 
can be kept small in comparison to its Roche-lobe, thus the tidal 
forces that circularise the system are weaker. By using the mech¬ 
anism of [Tout & Eggletoli| ( 1988| l, the wind mass loss depends 
on the orbital phase, with it being stronger at periastron than 
apastron. This difference between the mass lost at periastron and 
apastron can increase the orbital eccentricity, depending on the 
mass ratio of the system and the fraction of the wind mass loss 
that is accreted by the companion. 

The enhanced-wind-mass-loss model depends on only one 
parameter, Bwind, which Tout & Eggleton (19881 estimate at 


Bwind = 10"^ to explain the mass inversion in the pre-RLOE sys¬ 
tem Z-Her. The wind-mass-loss rate in function of the orbital 
phase 6 is: 


Mwind(^) ~ MReimers ' 


1 -H Bwind ■ min 




( 1 ) 


where R is the stellar radius and Rl{0) is the Roche-lobe radius at 
phase 6. A detailed description of the wind mass loss and accre¬ 
tion fractions in MESA is given in appendix |A.3| The strength 
of the eccentricity pumping will depend on the mass-loss rate, 
and the fraction that is accreted by the companion. We follow 
the model proposed by Soker ( 2000| l to calculate the change in 
eccentricity: 


^ml “ I 
Jg 


-H 2|Macc(6')| - - 

Md-tMa Ma 


(e + cos 0) do, (2) 


where Md and Ma are the donor and accretor mass. Moo is the 
mass lost at phase 0 to infinity, Macc is the mass lost at phase 0 
accreted by the companion and e is the eccentricity. Mass lost to 
infinity increases the eccentricity. Mass that is accreted, drives 
a change in eccentricity that is positive only if the donor mass 
is lower than the accretor mass. See appendix |A.4| for a detailed 
description of this mechanism. 

In Eig. the eccentricity as a function of the mass of the 
donor star is plotted for a binary model with initial period of 
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Fig. 3. The eccentricity in function of donor mass for models with dif¬ 
ferent values for wind-enhancement parameter fiwind: 100 (solid blue), 
5 000 (dashed red) and 10000 (dot-dashed green). The initial period is 
600 days, and the initial donor and companion masses are 1 . 2 - 1 - 0.8 Mq. 
For discussion see section[5] 


600 days, a donor mass and companion mass of 1.2 and 0.8 Mq, 
and different values of Bwind: 100, 5 000 and 10000. As can be 
seen, this method can result in a significant final eccentricity of 
the orbit after the RGB evolution of the primary, if the enhance¬ 
ment parameter is sufficiently large. For the system shown in 
Fig.[^ the circularisation is overcome when Bwind ^ 5000, while 
an sdB star is only formed when 10 < Bwind $ 100. The max¬ 
imum wind-mass-loss rates in these systems are: Mwind.max = 
10“^'^, 10“®'^, 10“® ° Mq yr“' for the systems with respectively 
Bwind = 100, 5 000 and 10000. For comparison, the wind-mass- 
loss rate for the same system without enhancement would be 

Mwind.max = 10“®'^ Mq yr“'. 

When comparing all calculated models, we find that the 
amount of mass that the donor needs to lose to maintain an ec¬ 
centric orbit is so high that the final core mass is too low for He 
ignition. The donor star will then end its life as a cooling He-WD 
on an eccentric orbit. To reach a final donor mass high enough 
to ignite He, the mass loss has to be lower so that the system 
completely circularises on the RGB. Therefore, the sdB binaries 
formed in this channel are all circularised. The final mass of the 
donor changes with the initial parameters of the system, like or¬ 
bital period, donor and accretor mass, and wind accretion frac¬ 
tion, but none of the combinations of initial binary parameters 
resulted in an eccentric system containing an sdB companion. 
Siess et aT|(|2014 1 gives a good overview of the effect of the ini¬ 


tial binary parameters and the wind-enhancement parameter on 
the final parameters of the binary system. 

We conclude that in the parameter regime considered here, 
the progenitors of sdBs do not evolve with a strong enhanced 
mass loss on the RGB. To form an sdB by tidally-enhanced wind 
mass-loss, the enhancement of the wind needs to be small, so 
essentially no eccentricity pumping occurs. Furthermore, if there 
is a significant enhanced wind mass-loss, the sdB progenitor will 
lose its hydrogen envelope in a stellar wind, and will not undergo 
RLOF during its later evolution. 


6. phase-dependent RLOF 

When there is no strong enhanced wind mass loss during the 
RGB, the sdB progenitor will eventually fill its Roche-lobe and 
start RLOF. If this mass loss happens on a slightly eccentric or¬ 


bit, the mass-loss rate will not be constant over the orbit, and the 
mass loss can have an eccentricity-pumping effect. The strength 
of the eccentricity pumping will depend on the mass-loss rate, 
and the fraction that is accreted by the companion, in the same 
way as in the tidally-enhanced-mass-loss mechanism. 


6.1. Model and input parameters 


In this section we will describe behaviour of the mass lost from 
the donor star in our study of phase-dependent RLOF. This mass 
can be accreted by the companion star, or is lost from the system 
to infinity. To describe the mass lost from the system to infinity, 
the method of |Tauris & van den Heuvel| ( |2006| l is used. A de¬ 
tailed description of that mechanism is given in appendix |A.2| 
The total mass loss from the system is subdivided in three frac¬ 
tions: from around the donor star (a), which carries the angu¬ 
lar momentum of the donor; mass transferred to the vicinity of 
the companion through the inner Lagrange point, and lost from 
around the companion as a fast wind (J3). This lost mass carries 
the angular momentum of the companion star; And through the 
outer Lagrange point (6) which is modelled as a circumbinary 
toroid with radius 1.25 times the binary separation (Rtoroid = 
y - 1.12, based on |Pennington| 1 9^ . The actual fraction of the 
mass loss accreted by the companion is defined by the aforemen¬ 
tioned fractions as:e = I - a - /3 - 6. 

The default system has an initial period of 700 days, initial 
eccentricity of 0.3, an sdB progenitor mass of 1.2 Mq and a com¬ 
panion mass of 1.0 Mq. The companion star will be on the main 
sequence during the relevant part of the evolution, thus the ra¬ 
diative dissipation mechanism for the tidal energy is assumed. 
The sdB progenitor will be on the RGB when the tidal forces are 
strongest, and the dissipation mechanism is that of a convective 
star. Dissipation becomes radiative when the envelope is com¬ 
pletely lost. A maximum value for the RLOF mass-loss rate is 
set at 10“^ Mq yr“\ the mass-loss fractions - jS- 0.35, 6 - 
0.30 and the location of the outer Lagrange point is Rl 2 = 1 -25 a, 
thus y - 1.12. The accretion fraction onto the companion star is 
zero in this default model. These parameters are summarised in 
Table Only the parameters mentioned in the text are changed. 

As it is not known when exactly a common envelope would 
start to form we assume that RLOF is stable in our models. To 
this end, we have capped the mass-loss rate during RLOF at 
10“^ Mq yr"' and applied an ad-hoc upper limit on the Roche- 
lobe overfilling of max(R/RL) = 1.25, meaning that models in 
which the donor star radius exceeds 1.25 ■ Rl are discarded. The 
latter limit is also imposed as our physical model breaks down 
at high Roche-lobe overfilling. It is still possible that in these 
circumstances a CE would form during RLOF. However, |Nele-| 
mans et al. ( 2000) l argues that if a CE would form during unstable 
RLOE in systems with a long orbital period, the CE would co¬ 
rotate with the binary, and there would be no friction between 
the binary and the CE, and consequently no spiral-in phase. Eur- 
thermore, based on BPS studies, [Clausen et al.| ( |2012| ) finds that 
the existence of long-period sdB systems with MS components 
indicates that RLOE at the tip of the RGB is stable, or otherwise, 
that a CE phase without spiral-in does exist. 


6.2. Parameter study 
6.2.1. Eccentricity evolution 

The evolution of the default model is plotted in Eig. [^ There are 
three events indicated on the figure, a: e^i > etidai, b: e^i < Stidai 
and c: log(etidai) < -15. Where is the eccentricity pumping of 


Article number, page 5 of 


18 





















A&A proofs: manuscript no. eccentric_sdbs_mesa 


Table 2. Standard parameters of the binary models. For the models fo¬ 
cusing only on phase-dependent RLOF, the disk mass is set to 0. The 
symbol column refers to the symbols used in the figures. 


Parameter 

Symbol 

Value 

sdB progenitor mass (Mq) 

Md 

1.2 

companion mass (Mq) 

Ma 

1.0 

period (d) 

p 

700 

eccentricity 

e 

0.3 

minimum eccentricity 

/ 

0.001 

Tidal forces 

dissipation type sdB 

/ 

Convective 

dissipation type companion 

/ 

Radiative 

Mass loss 

maximum M (Mq yr“') 

/ 

10^2 

fraction lost around companion 

a 

0.35 

fraction lost around donor 

p 

0.35 

fraction lost through L2 

5 

0.30 

location of L2 

7 

1.12 

accreted fraction" 

e 

0 

CB disk 

maximum mass (Mq) 

^disk 

0.01 

life time (yr) 

T 

10^ 

viscosity 

an 

0.01 

mass distribution 

cr(r) 

r-i 


Notes. The accreted fraction is defined by the other mass loss frac¬ 
tions as: e = I - a -- 6. 


the mass loss in phase-dependent RLOF and eudai is the circu- 
larisation due to stellar tides. These events indicate two phases 
in the period-eccentricity evolution of the binary. The timescale 
between these phases differs on the figure, and within one phase 
the time scale is linear. Underneath the figure the duration of 
each phase is given. 

Phase a-b: The sdB progenitor is close to completely filling its 
Roche-lobe (R/Rl = 0.97), and the mass loss rate reached 
10“^ ^Mq yr“', at which point the eccentricity-pumping 
forces are more effective than the tidal forces. Due to the 
increasing mass-loss rate, the eccentricity pumping stays 
stronger than the tidal forces, even though the latter increase 

due to further Roche-lobe overfilling (max R/Rl = 1.09). The 

eccentricity starts increasing, and reaches a value of 0.061 
by the end of this phase, an increase by a factor 60 compared 
to the assumed minimum eccentricity. During the mass-loss 
phase the period increases as well, from 847 to 945 days. 
When the star starts shrinking again, the mass-loss rate de¬ 
creases as well as the Roche-lobe-filling factor, until event 
‘b’, when the tidal forces again become higher than the ec- 
centi'icity pumping (M = 10“^ * Mq yr“*, R/R^ = 0.85). This 
phase of strong mass loss only takes 415 years. 

Phase b-c: When the RLOF diminishes, the tidal forces take 
over, and the binary starts circularising again for another 
16700 years. During this phase the sdB progenitor is con¬ 
tracting, thus the tidal forces are weakening. The eccentric¬ 
ity diminishes from 0.061 to 0.051, while the period slightly 
increases to 950 days. For this system, the eccentricity in¬ 
creased by roughly a factor 50. The evolution is plotted to 
the point where eudai < 10 after which there is no more 
significant change in eccentricity or period. 



415 yr 16700 yr 

Time 


Fig. 4. Time evolution of several binary properties during the RLOF 
phase. Three different events are indicated on the X-axes, a: e^ai > fitidai, 
b: e^i < fitidai and c: log(etidai) < -15. The time scale differs between the 
phases, but is linear withing each phase. The duration of each phase 
is shown under the figure. Panel A: the change in angular momentum. 
Panel B: mass loss through RLOF. Panel C: the tidal forces (blue) and 
eccentricity pumping due to mass loss (red). Panel D: the orbital eccen¬ 
tricity. Panel E: the orbital period. See section [6.2. l| for discussion. 

6.2.2. Effect of initial period and companion mass 

The effect of the initial period and companion mass is plotted in 
Fig.|6}A. With increasing initial period, the final period of the 
sdB binary will increase as well, while at the same time the ec¬ 
centricity decreases. A similar effect is visible when increasing 
the companion mass. The closer the companion mass is to the 
donor mass of 1.2 Mq, the larger the final orbital period, and the 
lower the eccentricity. We first discuss the effect of the initial bi¬ 
nary parameters, as they clearly illustrate the effect of the orbital 
period on the eccentricity pumping force that is also essential in 
the discussion of the process dependent parameters. 

To explain these effects, the evolution of several parameters 
for models with different companion masses of 0.8, 1.0 and 1.15 
Mq are shown in Fig.[^ These parameters are plotted in function 
of the donor mass instead of time, so that the different models 
can be more easily compared. The connection between eccen¬ 
tricity and orbital period is found in the change in mass-loss rate 
during RLOF. If the period of the system increases, the Roche- 
lobe overfilling is lower, varying from R/Rl <1.16 for the M^ = 
0.8 Mq system to R/Rl < 1.09 and 1.07 for M^ = 1.0 and 1.15 
Mq respectively. This lower overfilling of the Roche lobe leads to 
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a lower mass-loss rate and a shorter total time during which mass 
is lost at the maximum rate. For the models described here, the 
maximum mass-loss rate is the same, but the time during which 
it is sustained differs from roughly 25 years for the 0.8 Mq com¬ 
panion to 20 years for the 1.0 Mq companion and 15 years for the 
1.15 Mq companion. The difference in time scale seems small, 
but is significant when the mass-loss rate is 10“^ Mq yr“*. This 
diminishes the total eccentricity pumping force, even though the 
time during which it overpowers the tidal forces is longer for the 
1.15 Mq model (~ 380 yr) than for the 0.8 Mq model (~ 200 yr). 
Due to the lower Roche-lobe overfilling, the tidal forces will also 
decrease with increasing orbital period, with a factor 10 differ¬ 
ence between the 0.8 Mq model and the 1.15 Mq model. Even 
so, there is a net effect of lower eccentricity enhancement for the 
model with the heaviest companion. 

In Fig. only the period and mass range of models that re¬ 
sulted in an sdB binary are shown. For a system with a specified 
donor and accretor mass, the effect of the initial period on the fi¬ 
nal mass is important. If the orbital period is too short, the mass 
loss will be too strong, and the donor star will lose too much 
mass to ignite helium, ending up on the He-WD cooling track. 
The lower limit on helium ignition in MESA is around 0.45 Mq, 
slightly depending on composition. If the initial orbital period is 
too high, the mass-loss rate during RLOF will be too low, and 
the donor will ignite helium on the RGB and hence no sdB is 
created. In Fig. |^the HR diagram of four models with a donor 
and companion mass of 1.2 H- 1.0 Mq and initial periods of 600, 
650, 750 and 800 days are shown. The system with the lowest 
orbital period (panel A) ends up with a donor star mass of 0.448 
Mq, this is just under the He-ignition limit, and the donor star 
ends its evolution as a He white dwarf. The system with the 650 
days initial period (panel B) ends with a donor star of 0.456 Mq 
and ignites helium on the He-WD cooling track, which is called 
a late hot flasher. With an even higher period of 750 days, the fi¬ 
nal donor mass is 0.465 Mq, and the donor ignites helium shortly 
after departing of the RGB, during the evolution at constant lu¬ 
minosity, an early hot flasher. If the period is increased further, 
to 800 days, the donor will be too massive, and will ignite he¬ 
lium on the RGB when it still has a total mass of 0.929 Mq, 
after which mass loss is stopped until the donor finishes core-He 
burning and enters the AGB phase. 

6.2.3. Effect of mass-loss fractions and accretion 

The different ways to lose mass to infinity will change the final 
orbital parameters of the system. These parameters mainly in¬ 
fluence the amount of angular momentum that is removed from 
the system with the lost mass. The effect of these mass-loss frac¬ 
tions and the location of the outer Lagrange point is shown in 
figures I^B andj^C. In Fig. |^B, models for two different com¬ 
panion masses 1.0 and 1.15 Mq at a constant 6 - 0.2 are plot¬ 
ted, while changing the fractions of mass lost around the donor 
(a) and companion {jS). There is no accretion in these models 
{a + p + 6 - 1). When most mass is lost from around the donor 
the resulting period will be lower than when the mass is lost from 
around the companion star. This is easily explained by eq. |A.21| 
The mass that is lost from around the donor star will thus carry 
more angular momentum with it, which will result in a shorter 
orbital period. This change in period will influence the change in 
eccentricity by altering the mass-loss rates during RLOF, similar 
as was explained in Sect |6.2.2| From a certain threshold period 
that depends on the initial and mass-loss parameters, the eccen¬ 
tricity pumping is smaller than the tidal forces, and the orbits 
stay circularised. The effect of these two mass-loss parameters 



1.0 0.9 0.8 0.7 0.6 0.5 0.4 

Md (Mq) 


Fig. 5. Comparison of three models with different values for the com¬ 
panion mass: 0.8 (solid blue), 1.0 Mq (dashed red) and 1.15 Mq 
( dot-dashed green). Panel A: the Roche-lobe overfilling factor. Panel B: 
the mass loss rate during RLOF in Mq yr“'. Panel C: the tidal forces 
in log(s“*). Panel D: the eccentricity pumping forces in log(s“*). Panel 
E: net change in eccentricity per second. Panel F: mass ratio. Panel G: 
the eccentricity. Panel H: the orbital period in days. See section]^ for 
discussion. 

on the period is large, on the order of several hundred days. By 


Article number, page 7 of 


18 






















A&A proofs: manuscript no. eccentric_sdbs_mesa 




400 600 800 1000 1200 1400 

Period (days) 



Fig. 6. The effect of the initial and mass loss parameters of the RLOF models on the final period and eccentricity of the orbit. Panel A: Models 
with different initial periods (indicated by the black arrows) and different companion masses: Ma = 0.8 (blue stars), Ma = 1.0 (red dots) and Ma 
= 1.15 (green squares). Panel B: Models with different mass loss fractions a and /B, for two companion masses: 1.0 Mq (red circles) and 1.15 Mq 
( blue squares). The black arrow indicates decreasing a and increasing /3, S is constant at 0.2. The sum of the mass loss fractions is unity for each 
model. Panel C: Models with different mass loss fractions S: 0.1, 0.2 and 0.3, with y = 1.12 (red squares) and different values for y: 0.9, 1.0, 1.2 
and 1.3 with S = 0.3 (blue dots). These models don’t have accretion, thus a = yS = (1 - S)/2. Panel D: The effect of the accreted fraction for two 
companion masses and accretion fractions of 0%, 10% and 20%. = 1.0 Mq (red circles) and Ma = 1.15 Mq (blue squares). See Sects. [5^2^and 

l6.2.3l for discussion. 


changing from most mass lost around the donor to most mass 
lost around the companion, the final period can double. 


By increasing the fraction of mass that is lost through the 
outer Lagrange point (i5), more angular momentum is lost than 
if that mass was lost from around either of the binary compo¬ 
nents. The decrease in final period with increasing i5 then fol¬ 
lows directly from eq. A.21| Increasing the size of the circumbi- 
nary toroid representing the location of the outer Lagrange point 


(y) has the same effect, for exactly the same reason. Fig. |^C 
shows the effect of changing d and y. The models with vary¬ 
ing d have no accretion, thus a - j3 - (1 -6)12. The increase 
of the eccentricity with increasing 6 and y is directly related to 
the decrease of the orbital period, and thus a higher mass-loss 
rate. The two models on Fig. [^C with y values of 1.2 and 1.3 
result in binaries with very high Roche-lobe-overfilling factors 
{RjRh > L6). Note that our description breaks down at high 
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Roche-lobe-overfilling values, and the models in Fig.|^C with 
high values for y suffer from extrapolation. The effect of y is 
only given as an indication, as we keep it at 1.12 to represent the 
expected location of the outer Lagrange point. By changing 6 by 
only a small amount, the final period can again change drasti¬ 
cally, on the order of hundreds of days. Similar to the a and fi 
fractions, there is a certain threshold for 6 under which the orbit 
will stay circularised. The exact threshold value depends on the 
other parameters. 

The effect of the mass-loss fraction accreted onto the com¬ 
panion (e) is shown in panel D of Fig. By increasing the ac¬ 
creted fraction, the eccentricity increases while the orbital pe¬ 
riod stays more or less constant. By accreting a certain frac¬ 
tion of the mass loss, the size of the Roche-lobes will differ 
between the models. Higher accretion leads to slightly higher 
Roche-lobe overfilling, and a slightly higher eccentricity pump¬ 
ing. The period during which the eccentricity-pumping forces 
are stronger than the tidal forces also increases with increasing 
accretion rates. By increasing the accretion, the final eccentricity 
can almost be doubled, while the period remains constant. 


7. Circumbinary Disks 

The models with phase-dependent RLOF can indeed explain a 
certain part of the period-eccentricity diagram, but have prob¬ 
lems in the high-period high-eccentricity range, and cannot re¬ 
produce the circular systems. Circumbinary disks (CB disks) 
could potentially explain the high-period, high-eccentricity sys¬ 
tems as they add extra eccentricity-pumping forces on top of 
those from phase-dependent RLOF. 


7. 1. Model and input parameters 

CB disks can form around binaries during the Roche-lobe over¬ 
flow phase, if part of the mass can leave the system through 
the outer Lagrange points and form a Keplerian disk around 
the binary. The CB disk-binary resonant and non-resonant in¬ 
teractions have been described by Goldreich & Tremaine \\919\ 


and |Artymowicz & Lubow| ( 1994} , by using a linear-perturbation 
theory. The effect of the CB disk-binary resonances on the or¬ 
bital parameters has been the subject of many studies. In MESA 
we follow the approach outlined by Artymowicz & Lubow| 
P994| l and |Lubow & Artymowicz| ( |1996 l. 

The effect of the disk on the binary separation is given by: 


a 21 Jd 1 

a m Jb Ty 


(3) 


where To is the angular momentum of the disk and Jb the orbital 
angular momentum of the binary. / and m are integer numbers in¬ 
dicating the resonance with the strongest contribution a nd t„ i s 


the viscous evolution timescale (where Ty ~ see eq. A.45 i. 
The change in eccentricity depends on the change in binary sep¬ 
aration as: 


^disk “ ^ 

6 + - 

lOOe 


\m 


1 \ a 

a’ 


(4) 


for small eccentricities (e < 0.2), and decreases with 1/e for 
higher eccentricities. The derivation and implementation of these 
equations is given in appendix |A.6| 

The CB disk is described by its maximum mass (Mdisk), the 
mass distribution, the inner and outer radius, relative thickness 
near the inner rim (HIR), the viscosity of the matter (ao) and the 


total life time of the disk (t). In our model, the disk is formed 
by matter lost from the binary through the outer Lagrange point, 
and the disk itself loses mass at a rate determined by its life time. 
The mass in the disk is not constant, and only the maximum disk 


mass is a defined input parameter. The model of Artymowicz & 
Lubow| ( 1994i is only valid for a thin disk. Thus, H/R = 0.1 is 


assumed. The inner radius of the disk is determined based on 
smooth-particle-hydrodynamic (SPH) simulations, and the dust- 
condensation radius of the binary (see eq. |A.52 


A.54i. Based 


on observations and the assumed surface density behaviour, the 
outer radius of the disk is fixed at 250 AU (see also appendix 
lAAl. 


Maximum mass, life time, viscosity and distribution, are in¬ 
put parameters in the model. Based on observations of post-AGB 
disks dGielen et al.||20TT] [Hillen et al.||20T4l |Hillen|[M5l ), fol¬ 
lowing parameter ranges are assumed. The total disk mass can 
vary between 10^“* and 10“^ Mg. Disk life times after the end of 
RLOF range from 10"^ to 10^ years. The viscosity ranges from 
ao - 0.001 - 0.1, with the more likely range being oro < 0.01. 
The surface mass distribution decreases with increasing distance 
from the centre: cr(r) ~ r®, where -2 < D < -1. In the default 
model given in Table]^ a maximum disk mass of 10“^ Mg, a life 
time of 10^ years, a viscosity of 0.01 and a surface distribution 
of cr{r) ~ r~^ are chosen. 


7.2. Parameter study 

7.2.1. Eccentricity evolution 

In Fig.[7]the evolution of the CB disk mass, the period and ec¬ 
centricity, the change in angular momentum and the change in 
eccentricity are plotted in function of time. The model that is 
shown has the same parameters as the model displayed in Fig.|^ 
with the addition of a CB disk with a life time of 10^ yr. These 
parameters are also given in Table]^ Four different events are in¬ 
dicated on the figure, a. ^Jisk ^ ^tidaU ^ml ^ ^tidaU t:. ^ml ^ ^tidal 
and d: ejisk < etidai- These events define three different phases in 
the period-eccentricity evolution of the binary. On the figure, the 
time of each phase is linear, but the timescales between phases 
differ. The duration of each phase is plotted under the figure. 

Phase a-b: Mass lost from atmospheric RLOF is filling the CB 
disk, and the disk - binary interactions are strong enough 
to overcome the tidal forces, thus eccentricity starts to in¬ 
crease. The mass in the disk continues to grow while the 
donor star continues to expand, eventually filling and over¬ 
filling its Roche-lobe. The tidal forces continue to increase 
as well, and when the CB disk reaches its maximum mass 
after roughly 6100 years, the tidal forces again overtake the 
disk-binary interaction. By this time the eccentricity of the 
system reached 0.016 and the period decreased from 878 to 
852 days. The change in eccentricity is solely due to the disk¬ 
binary interactions, but the loss of angular momentum from 
the binary is caused mainly by the mass loss. During the re¬ 
maining ~1000 years when edi.sk < etidab there is a negligible 
amount of circularisation. 

Phase b-c: The disk reached its maximum mass, and is main¬ 
tained by RLOF. Due to phase-dependent mass loss, the ec¬ 
centricity continues to increase. What happens in this phase 
is very similar to what is shown in Figj^ However, because 
the binary has a higher eccentricity than that in the model 
without a disk, the eccentricity pumping due to mass loss is 
stronger. This leads to an eccentricity of 0.146, a tenfold in¬ 
crease. The orbital period first decreases and then increases 
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Fig. 7. Time evolution of several binary properties during the life time 
of the CB disk. Four different events are indicated on the x-axes. a: fidisk 
> e,idai, b: e„i > e,idai, c: e„i < e,idai and d: gdisk < eudai- The time scale 
differs between the phases, but is linear within each phase. The duration 
of each phase is shown under the figure. Panel A: fhe mass in fhe CB 
disk as a percentage of the maximum CB disk mass (0.01 Mq). Panel 
B: the change in angular momentum due to mass loss (red dashed line) 
and the CB disk - binary interaction (green dashed dotted line). Panel C: 
the tidal forces (blue full line), eccentricity pumping due to mass loss 
(red dashed line) and eccentricity pumping through CB disk - binary 
interactions (green dashed dotted line). Panel D: the eccentricity. Panel 
E: the orbital period. See section|7]for discussion. 


again, reaching 880 days. This whole phase lasts only 310 
years. 

Phase c-d: RLOF has ended, and the mass in the CB disk starts 
to decrease linearly. The CB disk-binary interactions con¬ 
tinue to increase the eccentricity, while at the same time an¬ 
gular momentum is transported from the binary to the CB 
disk, thus decreasing the orbital period again. The effect of 
the CB disk-binary interactions diminishes due to a decrease 
in disk mass while simultaneously, the resonances that are re¬ 
sponsible for the eccentricity pumping become less effective 
at higher eccentricities. When the disk is completely dissi¬ 
pated, the binary has an eccentricity of 0.245 and a period of 
855 days. 

We note that the transport of angular momentum from the binary 
to the disk will continue to take place as long as the disk lives, 
contrary to the change in eccentricity that has an upper limit of e 
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~ 0.7. However, no models we have tried reach this eccentricity 
limit. 


7.2.2. Effect of disk properties 


There are four parameters in the CB disk model that can influ¬ 
ence the CB disk-binary interactions: the maximum mass in the 
disk, the life time of the disk, the viscosity parameter and the as¬ 
sumed distribution of the mass with radius. Next to the disk pa¬ 
rameters, the rate at which mass is fed to the disk is determined 
by mass-loss fraction 6. The effect of these five parameters is 
shown in Fig.[^ panels A, B, C and D. For each of maximum 
disk mass, viscosity, mass distribution and 6, models with five 
different disk life times are shown. These life times can also be 
interpreted as the evolution in the period-eccentricity diagram af¬ 
ter the end of RLOF. Based on the equations that govern the CB 
disk-binary interaction given in section [A!6] the effect of each pa¬ 
rameter can be easily explained. When the disk survives longer, 
the interaction with the binary will last longer, and the increase 
in eccentricity / decrease in period will be stronger. 

By increasing the maximum disk mass as shown in panel A, 
the final eccentricity will be higher, while the final period will be 
slightly lower. The mass-loss rates during RLOF are high enough 
that there is very little difference in time to fill a disk of 0.001 
Mq or 0.01 Mg. An increased disk mass will lead to a higher 
angular momentum of the disk, and, according to eq. A.44 to a 
higher change in binary separation and eccentricity. 

The viscosity of the disk (cro) shown in panel B, also has a 
linear effect on the change in binary separation and eccentric¬ 
ity (see eq. |A.45[ ). A higher viscosity will lead to a higher final 
eccentricity and lower period. 

The effect of the mass distribution in the disk is plotted in 
panel C. The mass distribution will determine the mass fraction 
that resides close to the inner rim of the disk, where it has the 
strongest effect on the binary. Normally a radial mass distribu¬ 
tion of r ' is chosen, which takes into account that the disk is not 
just flat, but also has a specific thickness that varies with the ra¬ 
dius. A distribution of assumes that the thickness of the disk 
is more constant with radius. Under that assumption, there will 
be more mass closer to the inner rim, and the CB disk-binary in¬ 
teractions will be stronger, again leading to a higher eccentricity 
and lower period. 

The rate at which mass is fed to the CB-disk is defined by 
mass-loss parameter 6. The effect of changing this parameter is 
displayed in panel D. As the disk will reach its maximum mass 
earlier, the eccentricity-pumping effect of the disk will slightly 
increase with increasing 6. However, the main effect of this pa¬ 
rameter is a change in orbital period. As discussed in Sect. 6.2.3| 
increasing 6 will decrease the orbital period because more an¬ 
gular momentum is lost. This results in higher eccentricities 
reached during the RLOF phase. 


7.2.3. Effect of binary properties 

In panels E and F of Fig.[^ the effect of several binary param¬ 
eters is shown. These are also discussed in the Sect. 16.2.21 and 
are shown here, to indicate how strong their influence on the dis¬ 
tribution in the period-eccentricity plane can be. It is especially 
the range in period space that is altered by changing initial pa¬ 
rameters such as companion mass and initial period. Where the 
disk-binary interactions have an influence of roughly five per¬ 
cent on the final period, changing the initial binary parameters 
can result in a change of several 100 days. The main influence 
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Fig. 8. The effect of the CB disk and binary parameters on the final period and eccentricity of the orbit. All models are plotted for 5 different disk 
life times t = 1 • 10“*, 2.5 • 10"*, 5 • 10"*, 7.5 • 10"* and 1 • 10^ years. Panel A: Models with two different maximum disk masses: 10“^ Mq (blue squares) 
and 10“^ Mq (red circles). Panel B: Models with a different disk viscosities: = 0.1 (blue squares) and 0.01 (red circles). Panel C: Models with 

two different disk surface mass distributions: cr(r) ~ (blue squares) and cr(r) ~ (red circles). Panel D: Models with three different values 
for mass loss parameter 6: 0.1 (blue squares), 0.2 (red circles) and 0.3 (green stars). Panel E: Models with three different initial orbital periods: 
Pi = 600 d (blue squares), 650 d (red circles) and 700 d (green stars). Panel E: Models for three different companion masses: Ma = 0.8 Mq (red 
circles), 1.0 Mq (blue squares) and 1.15 Mq (green stars). See Sects. [7!T^ and [T2.3| for discussion. 

of the CB disk is then in reaching higher eccentricities. Where maximum 0.15, the models containing a CB disk can reach ec- 
the phase-dependent RLOF models could reach eccentricities of centricities of maximum 0.34. 
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Table 3. Parameter ranges used to determine the period-eccentricity dis¬ 
tribution of the phase-dependent RLOF models and the CB-disk mod¬ 
els. Not all parameter combinations result in an sdB binary. 


Parameter 

Range 

Binary 

sdB progenitor mass (Mq) 

1.0-1.50 

companion mass“ (Mq) 

0.8-1.45 

period (d) 

500 - 900 

Mass loss 

a 

0.0- 1.0 

p 

0.0- 1.0 

6 

0.0 - 0.30 

y 

1.12 

e 

0.0 - 0.25 

CB disk 

maximum mass (Mq) 

10-4 - 10-2 

lifetime (yr) 

104 - 10^ 

viscosity 

0.001 - 0.01 

mass distribution 

r-i 


Notes. The companion mass is always lower as the donor mass. 


8. Period-eccentricity distribution 

By varying the initial binary parameters described in Sect. |4] an d 
the model-dependent parameters described in Sects. and[^ a 
significant area of the period-eccentricity diagram can be cov¬ 
ered. In Fig. both the region covered by models with only 
phase-dependent RLOF (red shade with solid border) and that 
covered by models including a CB disk (green shade with dashed 
border) is plotted on top of the observed systems (blue circles). 
The parameter ranges used to obtain this distribution are given 
in Table]^ As the models with only tidally-enhanced wind mass- 
loss cannot produce eccentric sdB binaries, they are not shown 
here. We see that by including eccentricity-pumping processes 
during the evolution, we are able to form models with sdB pri¬ 
maries at wide eccentric orbits. 

From Fig.[^ we can conclude that the region covered by the 
models with phase-dependent RLOF does not correspond with 
the observations. The effect of most parameters is to reduce the 
orbital period, while increasing the eccentricity, leading to an ec¬ 
centricity pumping effect that is most efficient at short periods. 
At high orbital periods, the eccentricity pumping is too weak to 
overcome the circularisation, leading to circular orbits at peri¬ 
ods over ~1250 days. It is thus possible to recreate the short- 
period systems with moderate eccentricity, but the highly eccen¬ 
tric long-period systems cannot be reached. Another discrepancy 
is that the circular systems at orbital periods around 750 days 
cannot be reproduced either, when the eccentricity-pumping ef¬ 
fects are active. The only way to model these circular systems 
at low orbital periods is by turning off all eccentricity-pumping 
mechanisms. 

The models including a CB disk have a larger coverage than 
the models with only phase-dependent RLOF. However, even 
the disk models cannot completely reproduce the observed sys¬ 
tems. They suffer the same problems as the models with only 
phase-dependent RLOF. Again, the effect of almost all param¬ 
eters is to increase the eccentricity at the cost of lowering the 
orbital period. Even though higher periods and eccentricities can 
be reached, the system with the highest orbital period and eccen¬ 
tricity can still not be covered. 
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Fig. 9. The approximate region of the period-eccentricity diagram that 
can be explained by the RLOF models (red shaded region with solid 
border) and models containing a CB disk (green shaded region with 
dashed border). The observed systems are plotted in blue dots. 

The period-eccentricity distribution of the two eccentricity 
pumping models seems to indicate that there is a trend of higher 
eccentricities at lower periods, opposite to the observed trend of 
higher eccentricities at higher periods. This is under the assump¬ 
tion that the method-depending parameters might vary indepen¬ 
dent of each other and the initial binary parameters. However, it 
is not unreasonable that for example the mass-loss fractions (a, 
P and 5) depend on other binary properties like the mass ratio. 
Such a dependence between different parameters might indeed 
reproduce the observed period-eccentricity trend, but the models 
presented here can not prove nor disprove any such dependence. 

9. Summary and conclusions 

The goal of this article was to test if three eccentricity-pumping 
mechanisms proposed in the literature, could recreate the ob¬ 
served period-eccentricity distribution of long-period binaries 
containing a hot subdwarf B star. The three proposed mech¬ 
anisms are 1. tidally-enhanced wind mass-loss, 2. eccentricity 
pumping through phase-dependent RLOF and 3. the interac¬ 
tion between the binary and a circumbinary disk formed dur¬ 
ing RLOF combined with phase-dependent RLOF. The tidally- 
enhanced wind mass-loss method is not combined with any of 
the two other methods, as the point of the enhanced wind is to 
prevent the donor star from filling its Roche lobe, and in that way 
reducing the tidal forces. This would be counterproductive in the 
context of phase-dependent Roche-lobe overflow. 

In this work we attempted to describe the effects of model pa¬ 
rameters on the relationship between final period and eccentric¬ 
ity. Further work must be done in order to combine these mod¬ 
els with realistic distributions of starting parameters in a binary- 
population-synthesis context. 

Creating eccentric orbits with a tidally-enhanced-mass-loss 
mechanism is based on two processes. By increasing the wind 
mass loss, the radius of the donor star is kept well within its 
Roche lobe, thus reducing the tidal forces and maintaining the 
initial eccentricity of the orbit. The tidally-increased wind mass- 
loss depends on the orbital phase. This phase-dependent mass 


18 













J. Vos et al.: Testing eccentricity pumping mechanisms to model eccentric long period sdB binaries with MESA 


loss will further increase the eccentricity. The models calculated 
with MESA indeed lead to binaries on an eccentric orbit. How¬ 
ever, the amount of mass that the donor star needs to lose is so 
large that the final mass of the donor star is lower than the mini¬ 
mum core mass necessary to ignite helium. All eccentric binaries 
formed through this channel end up as He-WDs, while all sdB 
stars ended on a circularised orbit as the wind-enhancement pa¬ 
rameter needs to be very low. 

The phase-dependent RLOF mechanism is based on the 
same principle as the phase-dependent wind mass loss. By as¬ 
suming a minimal eccentricity, the mass-loss rate will be higher 
during periastron than during apastron, and the eccentricity may 
increase. This method can create eccentric sdB binaries over a 
reasonable period-eccentricity range. The maximum eccentric¬ 
ity tends to be around 0.15, while the period range for eccentric 
sdBs is 600 to 1200 days. sdB binaries with longer orbital pe¬ 
riod, reaching ~1600 days can be formed but are circularised 
within our parameter range. This mechanism can explain the ob¬ 
served systems with moderate eccentricities (0.05 < e < 0.15) 
on shorter orbital periods (P < 1100 d). Neither the systems with 
higher eccentricities, nor the circularised low-period systems can 
be formed this way. 

Adding CB disks to the models with phase-dependent RLOF 
does further increase the eccentricity of the produced bina¬ 
ries due to resonances between the binary and the CB disk. 
This binaty-CB disk interaction increases the eccentricity while 
slightly decreasing the orbital period. These models can reach 
eccentricities up to 0.35 when assuming reasonable values for 
the disk parameters, and could potentially reach far higher ec¬ 
centricities if more extreme values for maximum disk mass or 
surface mass distribution would be adopted. The disk model can 
reproduce most of the observed sdB binaries, except the circular 
short-period systems. 

When the models are compared to the observed sample, it 
is clear that phase-dependent RLOF in combination with CB 
disks can cover a significant part of the period-eccentricity di¬ 
agram. However, there remain conflicts. The circular systems at 
short periods can not be formed, as well at the system with the 
higher period and eccentricity. Part of the discrepancy between 
the observed and modelled period-eccentricity diagram can be 
explained by insufficiently accurate models, while it is also pos¬ 
sible that certain areas of the period-eccentricity diagram have 
not yet been populated due to a lack of observations. 

The observed trend of higher eccentricities at higher orbital 
periods does not follow directly from the method-depending pa¬ 
rameters in the presented models. In our future research we will 
investigate whether certain parameters are dependent on the ini¬ 
tial or current configuration of the binary. Such a dependence 
could potentially result in the observed period-eccentricity trend. 
However, the models presented here can neither prove nor dis¬ 
prove any such trend. 

All three tested eccentricity-pumping mechanism are derived 
under certain assumptions, such as isotropic mass loss for ec¬ 
centricity pumping through phase-dependent mass loss, or the 
approximation of a thin disk for the CB disk-binary interactions. 
Furthermore, values for several input parameters, for example 
the mass loss fractions (a, (3 and y), are badly constrained or cur¬ 
rently unknown, and are also likely to be functions of the system 
parameters such as mass ratio and companion mass rather than 
constants. While stable Keplerian discs are likely common in 
other evolved binary systems (see introduction), they are not yet 
documented around sdB wide binaries. We advocate here a spe¬ 
cific search for circumstellar matter around wide sdB binaries. 


In this article we have shown that we are able to form bi¬ 
nary models with an sdB primary at wide eccentric orbits by 
eccentricity pumping of phase-dependent RLOF and CB disks. 
Tidally-enhanced wind mass-loss is unlikely to contribute to the 
formation of eccentric sdBs. Small discrepancies between the 
observed systems and the theoretical models remain, and will 
need to be addressed in future work. 
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Appendix A. 1: Evolution of the orbital parameters 

The orbital angular momentum of a binary system is: 


J orb “ a 


Ga(l-e^) 

I Md + M, ' 


(A.l) 


Where a is the binary separation, and e is the eccentricity of the 
orbit. Mj and Ma are respectively the mass of the donor and the 
accretor, and G is the gravitational constant. The evolution of 
the binary separation is obtained from the time derivative of Eq. 
(lATl: 


d _ 2 '^orb 

n J orb 


M, 


^ Md-bMa ^ ^ ee 


Md Ma Md + Ma 1 


,2 ■ 


(A.2) 


In MESA, the mass loss and accretion rates, the change in or¬ 
bital angular momentum, and the change in eccentricity are com¬ 
puted, and then used to update the orbital separation. 

The change in mass of both components depends on the mass 
lost in stellar winds, the mass lost through Roche-lobe overflow 
(RLOE), and the fraction of mass that is accreted by the com¬ 
panion. Eor the donor star the net mass loss is: 


Md — Mrlof + Mwind.d ^wind.a ' M^ind.ai 

while for the accretor: 


Ma — Mrvind.a ^wind.d ' M^ind.d ^rlof ' Mrlof. 


(A.3) 


(A.4) 


Where e^ind is the fraction of the wind mass loss that is accreted 
by the companion and eriof is the fraction of the RLOE mass loss 
accreted by the companion. The amount of mass lost to inhnity 
is then: 

44oo — (1 “ eriof) ■ 4ifrlof "b (1 ewind.d) ' ^ifwind.d 

■b (1 “ ewind.a) ' Mwind.a- (A.5) 

The change in orbital angular momentum has three main 
components: mass lost to inhnity in stellar winds, mass lost to 
inhnity during Roche-lobe overhow and the change in angular 
momentum due to resonances between the binary and a poten¬ 
tial circumbinary (CB) disk: 


•^orb “ '^rlof "b ./wind "b ./disk- 


(A.6) 


In the wide binary models considered here, the strongest con¬ 
tributions to 7orb will come from mass loss during Roche-lobe 
overhow and the interaction between the binary and a CB disk. 

The evolution of the eccentricity is governed by the tidal in¬ 
teractions between both stars, the effect of phase-dependent mass 
loss/transfer which can pump or decrease the eccentricity based 
on the mass ratio and the interaction between the binary and a 
circumbinary disk: 


^ ~ tides "b ^ml "b ^disk- 


(A.7) 


The relation between the orbital period and the separation is: 


Appendix A: Binary physics in MESA 

In this appendix we will list the physics used in the binary evo¬ 
lution module of MESA. Eor the exact software implementation 
we refer the reader to the open-source code, and future instni- 
ment papers. 


P^2n- 


(G(Md + MJ 


(A.8) 


All physical mechanisms relevant to the terms given in equa¬ 
tions (|A.3|-[A77 ]i are explained in the following subsections. 
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Appendix A.2: Roche-lobe overflow 

The instantaneous Roche-lobe radius of the donor star (Rl.ci) at 
a specific moment in the orbit is approximated by adapting the 
fitting formula for the effective Roche-lobe radius of [Eggleton 
( |1983| l by allowing the binary separation to depend on the orbital 
phase (6); 


RlAS) = 


0.49^^' 


2/3 


1 H- ecos(6l) 0.60 H- ln(l -H 


(A.9) 


Where a is the separation, e the orbital eccentricity. The mass 
ratios q^ and q^ are defined as: 


Md 




= —■ 

Ma Mi 


(A. 10) 


The Roche-lobe radius of the accretor can be obtained by invert¬ 


ing the mass ratio in Eq. A.9 thus replacing qi with q^. 

Mass-loss is implemented in MESA following the prescrip¬ 
tion of |Ritter| ( |1988| l and |Kolb & Ritter| ( |1990| l. This formal¬ 
ism differs between mass-loss from the optically thin region 
(Ri < *^hat from the optically thick region (Rj > Ru.d)- 

Mass-loss from the former region is given by: 


An approximation of the distance between the centre of mass of 
the donor and the LI point in units of the binary separation (xl) 
is given by |Erank et al. ( 2002| l: 


3Cl = 0.5 - 0.227 logio(q'). 


(A. 17) 


In the case that the stellar radius significantly overfills its 
Roche-lobe, the Roche-lobe will lie in the optically thick region 
of the donor. The mass loss is then calculated via: 


A Irt 

Mthick = Mo + 27 rE(^)^ F,(r,).—dP. (A. 18 ) 

GMi ^ p 


Where Mq is given in Eq. (A.14i, Eli and Pph are respectively 
the pressure at LI and at the stellar photosphere. Both the tem¬ 
perature T and the mean molecular weight p depend on the pres¬ 
sure. Fj is a function that depends on the adiabatic exponent Ti, 
and is given by: 




(A.19) 


Afthin = Mq ■ exp 


/ Rd - ^L,d \ 

\ /7f,Li /■ 


(A. 11) 


Where Rj is the donor radius and FIp^i is the pressure scale 
height at the inner Lagrange point (LI), which can be linked to 
the pressure scale height at the photosphere of the donor Hppi 
as: 


Hp^pi 1 RphR^ 

7(q) yiq) Pph g Mi' 


(A. 12) 


The mass loss equations depend on the eccentric anomaly 
through the Roche-lobe radius. The time steps that MESA uses 
in the evolution are much larger than one orbital period, thus to 
obtain the average mass loss rate to be used in a time step, the 
above equations are integrated over the orbit. 

RLOE is not necessarily conservative. To allow for mass 
to be lost from the system the formalism of|Tauris & van den| 


Heuvel (2006i and Soberman et al. (1997 1 is used. This system 


describes three fractions (a,/?, S) to lose mass from the system: 


Here Rph is the pressure at the donor photosphere and pph is the 
photospheric density of the donor. y{q) is a factor that depends 
on the mass ratio as: 


0.954 -H 0.025 logio(^) - 0.038 Wwiq), 0.04 < ^ < 1. 


yiq) = ^ 


0.954 -H 0.039 logio(^) -1-0.114 \og\iiq), I <q<20. 

(A. 13) 

Mq is the mass loss rate when the donor star just fills its Roche- 
lobe: 


In /RTeffd 
Mo = — 

ye \ Pph.d 


4d 


GMi 


Pph.d E(^o)- 


(A. 14) 


- a: mass lost from the vicinity of the donor as a fast wind 
(Jeans mode). This is modelled as a spherically symmetric 
outflow from the donor star in the form of a fast wind. The 
mass lost in this way carries the specific angular momentum 
of the donor star. 

- yS: mass lost from the vicinity of the accretor as a fast wind 
(Isotropic re-emission). A flow in which matter is transported 
from the donor to the vicinity of the accretor, where it is 
ejected as a fast isotropic wind. Mass lost in this way carries 
the specific angular momentum of the accretor. 

- d: mass lost from a circumbinary coplanar toroid. The radius 
of the coplanar toroid is determined by y as Rtoroid = y^n. 

The accretion efficiency of RLOE is then given by: 


Pph.d and pph,d are the mean molecular weight and density at the 
donor photosphere, Tefj.d is the effective temperature of the donor 
and R is the ideal gas constant. F(q) is a function of the mass ra¬ 
tio depending on the Roche geometry, given by (see also Meyer 
& Meyer-Hof'meister|1983 1: 


G(q) 


q 


ylgHq)-{'i +q)g{q) v^^.d 


(A. 15) 


friof-l-a-yS-d. (A.20) 

When the mode of mass loss is known, the change in angular 
momentum can be calculated. Mass accreted by the companion 
will not change the total angular momentum of the system, thus 
only fractions a,/? and 6 will have an influence on J. The total 
effect of mass loss through Roche-lobe overflow on the change in 
angular momentum is given by|Tauris & van den Heuvel](|2006|l: 


where g(q) is: 



1 

(1 -xiA' 


(A. 16) 


J rlof “ 


a + Pq^ + dy{l + qp M^af 
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(A.21) 
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Appendix A.3: Wind mass ioss 

The amount of mass lost through stellar winds is determined by 
any of the wind loss prescriptions defined in the stellar part of 
MESA. The binary module offers the possibility to boost this 
mass loss using the Companion Reinforced Attrition Process 
(CRAP) mechanism ( |Tout & Eggleton]|1988[ ), also known as 
tidally enhanced mass loss. In this model it is assumed that tidal 
interactions or magnetic activity are responsible for the enhance¬ 
ment of the wind. The enhancement is expected to have a similar 
dependence on radius over Roche-lobe radius as the torque in a 
tidal friction model, and there is an expected saturation when 
co-rotation is reached (in their model when R - Rl/2). 


^wind “ ^Reimers lid" ^wind ' min 


R 1 


RlI ’26 


(A.22) 


The factor Bwind is estimated by Tout & Eggleton ( 1988| l to be of 
the order Bwind ~ 10"*, but can vary significantly depending on 
which system needs to be explained. 

Part of the mass lost due to stellar winds is accreted by the 
companion. In the case of fast winds (Vwind » Vorb) the accretion 
fraction is given by the Bondi-Hoyle mechanism ( [Hurley et al.| 
2002^ . Here we give the equations for the donor, those of the 
accretor can be obtained by switching the d and a subscripts: 


^BH.d 




GMd 


ttBU 


1 


\ ^wind,a 7 


2fl2 (1 H- 


The velocities are: 


(A.23) 






^orb 

V^. . 
wind.a 


G(Md -H Ma) 


V, 


wind,a 


■ 2^iv ■ 


a 

GM, 


R. 


(A.24) 

(A.25) 

(A.26) 


The wind velocity used here is set proportional to the escape ve¬ 
locity from the stellar surface. Based on ob served wind veloc ities 
in cool supergiants (v^ind = 5-35 km s“*) (Kucinskas 1998 i, Pw 
is taken to be 1/8. The free parameter q-bh is set to 3/2 based on 
[Boffin & Jorisseri] ( |1988| ). 

The angular-momentum loss due to wind-mass loss, assum¬ 
ing a spherical-symmetric wind is: 


J wind.d — llfwind.d 




Mi + Ma 


2n r -T 

— Vl - e^. 

P 


(A.27) 


Eor the change in total angular momentum of the binary due to 
stellar winds, the contribution of the accretor need to be added 
and ^wind — 'fwind.d "b 'fwind.a- 


Appendix A.4: phase-dependent mass ioss 

When mass loss, weather due to stellar winds, RLOE or other 
reasons, is not constant during the binary orbit, it will have an ef¬ 
fect on the eccentricity of the system. There are multiple causes 
of the periodicity of mass loss and accretion. Eor example, the 
mass loss might be caused by stellar pulsations, which can trans¬ 
fer more mass at the periastra that coincide with a maximum 
stellar radius. The methods implemented in MESA are phase- 
dependent wind-mass loss through tidal interactions, and phase- 
dependent RLOE on eccentric orbits. The latter will boost mass 
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loss near the periastron passage while during apastron the star is 
completely inside its Roche lobe. The effect of phase-dependent 
mass loss on the orbital eccentricity was studied by [Soker ( 2000[ l 
based on the theoretical work of [Eggleton] P006| Ch. 6.5). 

In calculating the effect on the eccentricity we have to dis¬ 
tinguish between mass lost to infinity, and mass that is accreted 
by the companion. Assuming isotropic mass loss, the change in 
eccentricity for mass lost from the system is given by: 




(A.28) 


Where 6 is the true anomaly, and Mao(0) is the mass lost from the 
system at a specific phase 6 during the orbit. Erom this equation 
one can see that constant mass loss will not change the eccen¬ 
tricity, if however the fraction of mass lost near periastron is sig¬ 
nificantly larger than during the remainder of the orbit, Eq. [A.28[ 
predicts a positive e. 

The effect on e of mass accreted by the companion, again 
under the assumption of isotropic mass-loss, is given by: 

eacc(0) = 2|Macc(e)l ( " 4“ I (^'29) 

\Mi Mj 


Where is the mass accreted by the companion at a spe¬ 

cific phase 6 during the orbit. As mass transfer is expected near 
periastron, the eccentricity will increase if Mi < Ma, which is 
required in the case of stable mass transfer. 

To obtain the total change in eccentricity in one orbit, Eq. 


A.28 and A.29 are integrated over the orbit: 


f'^ 


^ml “ I [ ^lost(^) ^acc (0)]de. 


(A.30) 


Appendix A.5: Tidai forces 

The gravitational forces acting on the components in binary sys¬ 
tems induce a deformation of their structure, and creates tidal 
bulges. These bulges are lagging behind, thus the gravitational 
attraction generates a torque on those bulges, which forces the 
synchronisation and circularisation of the stellar and orbital rota¬ 
tion. The orbital parameters change because the star’s rotational 
energy is dissipated into heat. There are two cases one has to 
consider, stars with convective envelopes and stars with radiative 
envelopes. In the former, the kinetic energy of tidally-induced 
large-scale currents is dissipated into heat by viscous friction of 
the convective environment. In radiative stars it is mainly radia¬ 
tive damping of gravity modes that functions as dissipation pro¬ 
cess. See [Zahn] ( [2005] l for a review on tidal dissipation. 

We use the formalism developed by [Hut| ( fI98I[ l to calculate 
the effect of tides on circularisation and synchronisation. This 
formalism, developed based on the weak-friction model, results 
in (where we again give only the equations relevant for the donor 
star, and those of the accretor are obtained by switching the d and 
a subscripts): 


^tides.d — 27 




yr 


(A.31) 
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for the circularisation, while the synchronisation is given by; Appendix A.6: Circumbinary disks 


n,.des,d-3|-| -(-j 


d n 


(1 - e2)6 




Q 




yr 


-1 


Here / 2-5 are polynomials in e; 


^ 15 9 45 H 5ft 

U{e^) = 1 + 
fsie'^) = 1 + 3 eY 


(A.32) 

(A.33) 
(A.34) 
(A.35) 
(A. 36) 
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rad,d 
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^tides “ ^tides.d T ^tides.i 


Circumbinary disks (CB disks) can form around binaries dur¬ 
ing the Roche-lobe overflow phase, if part of the mass can leave 
the system through the outer Lagrange points and form a Ke- 
plerian disk around the binary. The CB disk - binary resonant 
and non-resonant interactions have been described by |Goldre-| 
ich & Tremaine \\919\ and Artymowicz & Lubow ( 1994[ ), by 
using a linear-perturbation theory. There are two main assump¬ 
tions in these models; the first is that the disk is thin (0.01 < 
HIR <0.1, where H and R are respectively the thickness and 
the half-angular-momentum radius of the disk). The second as¬ 
sumption is that the nonaxisymmetric potential perturbations are 
small around the average binary potential. 

The effect of the disk - binary resonances on the orbital pa¬ 
rameters has been the subject of many studies. In MESA we fol¬ 
low the approach of|Artymowicz & Lubow ( 1994|l and Lubow & 


and Tg is the radius of giration. Furthermore kjT depends on the 
star being convective or radiative. In the case of a convective 
envelope, k/T is given by Hurley et al. (2002 eq. 30-33), based 
on |Rasio et al.| ( [r996| l: 


(A.37) 


Where Tconv is the eddy turnover time-scale (the timescale on 
which the largest convective cells turn over): 


(A.38) 


Artymowicz| ( | 19961 1. Ths model of Lubow & Artymowicz ( 1996| l 

for small and moderate eccentricities (e < 0.2) is based on the re¬ 
sult of smooth particle hydrodynamic (SPH) simulations, which 
show that for a disk in which the viscosity is independent of 
the density, the torque is independent of resonance strength and 
width. The variation in the orbital separation due to the inner and 
outer Lindblad resonances is ( |Lubow & Artymowicz|199^ : 

(a.44) 

a m 7b Ty 

Where I and m are integer numbers indicating which resonance 
has the strongest contribution. 7 d and 7 b are respectively the an¬ 
gular momentum in the disk and in the binary, and r,, is the vis¬ 
cous evolution timescale, which is the timescale on which matter 
diffuses through the disk under the effect of the viscous torque. 
It is given by; 


1 


And /conv is a numerical factor depending on the tidal pumping 
time-scale Ptid: 


ao 


(r) 


(A.45) 


(A.39) 


(A.40) 


Where ao is the viscosity parameter of the disk and Q. 1 , is the 
orbital angular frequency. With a disk viscosity of an = 0.1, the 
viscous timescale is typically on the order of 10^ yr. 

To determine the angular momentum of a Keplerian disk, one 
need to know the surface mass distribution cr in the disk. 7 d is 
then given by: 

^^out 

Jd - I r ■ cr ■ V dA = I r ■ cr ■ s/CMsIr 2 n r dr, (A.46) 


= r r-o--vdA = r 

Ja 


Menv and Penv are respectively the mass in the outer stellar con¬ 
vective zone, and the stellar radius at the base of that zone. L is 
the luminosity. 

In the case of a radiative envelope, the damping is caused by 
a range of oscillations that are driven by the tidal field. kjT is 
then ( Hurley et al.|2002' eq. 42-43 based on Zahn|1973 1977| l: 


where the Keplerian velocity of an element in the disk at distance 
r from the centre of mass, when neglecting the mass of the disk 
is: V - -\JG{Md + M^)lr. Pin and Pout are the inner and outer 
boundaries of the disk. We assume the surface distribution in the 
disk to depend on the radius to the power of a free parameter 6: 




(A.47) 


(A.41) 


Where £2 is a second order tidal coefficient which can be fitted 
with: 


(A.42) 


The distribution constant Dc depends on the total disk mass and 
6. With this distribution the disk angular momentum can be cal¬ 
culated; 

pRom 

Jd = 2n yjGiMi + M^) dr. (A.48) 

■JR\u 

Which depending on parameter 6 has the following solutions: 


The total effect on the change in orbital eccentricity is then 
obtained by combining the tidal forces on both components: 


Jd — ' 


(A.43) 




, 5 / 2 - 
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•In Pin), b = 5/2. 
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(A.49) 


As shown by the SPH simulations of |Lubow & Artymowicz 
( |1996| l, only the inner part of the disk plays a role in the disk 
- binary interactions. Thus in the previous equations the angular 
momentum is only calculated from the region between Ri„ and 
six times the binary separation (Rout - 6 a). 

To calculate the distribution constant in Eq. ( |A.47| l, one has 
to know the total mass in the disk (Mdi.sk). The total disk mass 
can then be related to the distribution constant as: 


Do = 


Mdisk 




rl-S 


dr 


(A.50) 


Which depending on 6 has two solutions: 


Do = 


1 (2-(5)Mdisk 

out m 


1 


Mdisk 


In In Rn 


inR,, 


6 + 2. 


6^2. 


(A.51) 


The main difference between the implementation here and 
that of |Dermine et al.| ( [2013[ l is that the disk-surface-distribution 
parameter 6 in our model is consistent throughout the calculation 
of the disk angular momentum and the distribution constant Do, 
while [Pennine et al.| ( |2013| l used 6 - 2.5 in the former, while 
(5=1 was used in the latter. 

We consider two processes to determine the inner radius of 
the CB disk. The interaction between the binary and the disk will 
clear the inner part of the disk. SPH simulations performed by 
Artymowicz Sc Lubow ( |1994| l show that the inner radius depends 
on the Reynolds number of the disk gas ()R = (HIRy^a^) and 
the eccentricity. A fitting formula of their results was derived by 
jPermine et al.| ( |2013| ): 


^i„.sPH = 1.7--log()RV^)AU. 

O 


(A.52) 


The second limiting factor on the inner radius is the dust conden¬ 
sation temperature. Based on Pullemond Sc Monnier (2010 Eq. 
1 - 12 ), and adding an offset for the binary separation we obtain: 


Td + R-u 


a(MuLi 


L^u OKMuW + MdLa) 

J 47r (Tbol (Md -H Ma)(Ld + Ta) 


(A.53) 


Where Ld and La are respectively the luminosity of the donor and 
accretor, and crboi is the Stefan-Boltzmann constant. Toond is the 
dust condensation temperature which we take at 1500 K. The last 
factor in Eq. A.53| is the average binary separation weighted by 
the luminosity of both components. The actual inner radius of the 
disk is the maximum of that determined from SPH simulations 
and the dust condensation radius: 


Riu — max [Rin.SPH , f^in.dust ] • (A.54) 

There are no direct observations of the outer radius of a CB 
disk around a post-AGB or post-RGB binary. Post-AGB and 
post-RGB disks are likely similar to protoplanetary disks, thus 
it can be assumed that the surface mass distribution of a CB disk 
does not follow the same behaviour along the whole radius. The 
inner part of the disk follows cr(r) ~ r^^,6 = 1 - 2 , while the 
outer part has an exponential decline in cr. Eor the CB disk¬ 
binary interaction only the inner part of the disk is important. 
The outer disk radius here then represents the radius where the 
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exponential drop in mass distribution starts, and the maximum 
disk mass represents the mass in this inner part of the disk. We 


assume an outer disk radius of 250 At/. See for example Bujarra- 
jbal et al.| ( |2007|[2()13||2015| l. Keep in mind that this outer radius 
is not the same as the six time the separation radius that is used 
in calculating the effective disk angular momentum. 

The change in eccentricity due to Lindblad resonances can 
be given in function of d/a, and depends on the eccentricity. 
Eor small eccentricities this can be calculated analytically. In the 
range e < 0.1 -y/oo the m = I resonance dominates, while in the 
region 0.1 -^/od < e < 0.2 the m - 2, I - I resonance is dom¬ 
inant. Eor small e the analytic form of e is (see also Lubow Sc 
|Artymowicz|200()|l: 
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^disk 


(TD 

1 ^ 


1 


vr 


e < 0.2. 


(A.55) 


Eor average eccentricities the above equation is extrapolated 
cording to an efficiency decreasing with 1 /e: 


ac- 


edisk = (—+C 5 ) ■0. 
\ e la 


2<e< 0.7. 


(A.56) 


Where the constants are C 4 = -1.3652 and C 5 = -1.9504 so 


that Eq. A.56 smoothly connects to Eq. A.55 while going to 0 
at e = 0.7. Eor larger eccentricities resonances that damp the 
eccentricity start to dominate, thus Cdisk = 0 in the range e > 0.7 
(e.g. Roedig et al.|201 ij l. 

The change in angular momentum due to disk - binary inter¬ 
actions is: 


4sk = Tocb ■ (A.57) 

Where d/a is given by Eq. ( |A.44| l. 

The mass feeding the CB disk is the fraction of the mass 
lost from the binary system through the outer Lagrange point 
during RLOE. The CB disk mass-loss rate is determined by the 
maximum mass in the disk Mdisk,max, and the life time of the disk 
Tdisk- Both are input parameters in our model. The rate in which 
the disk loses/gains mass is then: 

Mdisk = b IMriofI - . (A.58) 

^disk 

In this implementation the disk will continue existing for a pe¬ 
riod Tdisk after RLOE stops. 
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